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Abstract: In the simple mean-field 5/5" and SIR epidemic models, infection 
is transmitted from infectious to susceptible members of a finite population 
by independent p— coin tosses. Spatial variants of these models are proposed, 
in which finite populations of size N are situated at the sites of a lattice 
and infectious contacts are limited to individuals at neighboring sites. Scaling 
laws for both the mean-field and spatial models are given when the infection 
parameter p is such that the epidemics are critical. It is shown that in all 
cases there is a critical threshold for the numbers initially infected: below the 
threshold, the epidemic evolves in essentially the same manner as its branching 
envelope, but at the threshold evolves like a branching process with a size- 
dependent drift. 



1. Stochastic epidemic models 
1.1. Mean-field models 

The simplest and most thoroughly studied stochastic models of epidemics are mean- 
field models, in which all individuals of a finite population interact in the same 
manner. In these models, a contagious disease is transmitted among individuals of 
a homogeneous population of size N. In the simple SIS epidemic, individuals are 
at any time either infected or susceptible; infected individuals remain infected for 
one unit of time and then become susceptible. In the simple SIR epidemic (more 
commonly known as the Reed- Frost model), individuals are either infected, suscep- 
tible, or recovered; infected individuals remain infected for one unit of time, after 
which they recover and acquire permanent immunity from future infection. In both 
models, the mechanism by which infection occurs is random: At each time, for any 
pair (i, s) of an infected and a susceptible individual, the disease is transmitted 
from i to s with probability p = pn- These transmission events are mutually in- 
dependent. Thus, in both the SIR and the SIS model, the number Jt+x — of 
infected individuals at time t + 1 is given by 
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where St = is the number of susceptibles at time t and the random variables £ s 
are, conditional on the history of the epidemic to time t, independent, identically 
distributed Bernoulli-1 — (1 — p) Jt . In the SIR model, 



where Rt is the number of recovered individuals at time t, while in the SIS model, 



In either model, the epidemic ends at the first time T when Jt — 0. The most basic 
and interesting questions concerning these models have to do with the duration T 
and size ^2 t <T Jt °f the epidemic and their dependence on the infection parameter 
Pn and the initial conditions. 

1.2. Spatial SIR and SIS epidemics 

In the simple SIS and SIR epidemics, no allowance is made for geographic or social 
stratifications of the population, nor for variability in susceptibility or degree of 
contagiousness. Following are descriptions of simple stochastic models that incor- 
porate a geographic stratification of a population. We shall call these the (spatial) 
SIS—d and SIR—d epidemics, with d denoting the spatial dimension. 

Assume that at each lattice point x £ 7L d is a homogeneous population of N x in- 
dividuals, each of whom may at any time be either susceptible or infected, or (in the 
SIR variants) recovered. These populations may be thought of as "villages" . As in 
the mean-field models, infected individuals remain contagious for one unit of time, 
after which they recover with immunity from future infection (in the SIR variants) 
or once again become susceptible (in the SIS models). At each time t = 0, 1, 2, . . . , 
for each pair (i x , s y ) of an infected individual located at x and a susceptible indi- 
vidual at y, the disease spreads from i x to s y with probability a(x, y). 

The simple Reed- Frost and stochastic logistic epidemics described in section 11.11 
terminate with probability one, regardless of the value of the infection parameter 
p, because the population is finite. For the spatial SIS and SIR models this is 
no longer necessarily the case: If Ylx£Z d ^ x = 00 then, depending on the value 
of the parameter p and the dimension d, the epidemic may persist forever with 
positive probability. (For instance, if N x — 1 for all x and a(x, y) = p for nearest 
neighbor pairs x,y but a(x, y) = otherwise, then the SIS —d epidemic is just 
oriented percolation on Z d+1 , which is known to survive with positive probability 
if p exceeds a critical value p c < 1 [6J.) Obvious questions of interest center on how 
the epidemic spreads through space, and in cases where it eventually dies out, how 
far it spreads. 

The figure below shows a simulation of an SIS-1 epidemic with village size 
N x = 20224 and infection parameter 1/60672. At time there were 2048 infected 
individuals at site 0; all other individuals were healthy. The epidemic lasted 713 
generations (only the first 450 are shown). 

1.3. Epidemic models and random graphs 

All of the models described above have equivalent descriptions as structured ran- 
dom graphs, that is, percolation processes. Consider for definiteness the simple SIR 
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R t +i = Rt + Jt and 
Si+i = St — Jt+i, 
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St+i — St + Jt — Jt+i- 
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Fig 1. 



(Reed-Frost) epidemic. In this model, no individual may be infected more than once; 
furthermore, for any pair x, y of individuals, there will be at most one opportunity 
for infection to pass from x to y or from y to x during the course of the epidemic. 
Thus, one could simulate the epidemic by first tossing a p— coin for every pair x, y, 
drawing an edge between x and y for each coin toss resulting in a Head, and then 
using the resulting (Erdos-Renyi) random graph determined by these edges to de- 
termine the course of infection in the epidemic. In detail: If Yq is the set of infected 
individuals at time 0, then the set Y\ of individuals infected at time 1 consists of 
all x Yq that are connected to individuals in Yq, and for any subsequent time n, 
the set Y n +\ of individuals infected at time n+1 consists of all x ^ Uj< n Yj who are 
connected to individuals in Y n . Note that the set of individuals ultimately infected 
during the course of the epidemic is the union of those connected components of 
the random graph containing at least one vertex in Yq . 

Similar random graph descriptions may be given for the simple SIS and the 
spatial SIS and SIR epidemic models. 

1.4- Branching envelopes of epidemics 

For each of the stochastic epidemic models discussed above there is an associated 
branching process that serves, in a certain sense, as a "tangent" to the epidemic. 
We shall refer to this branching process as the branching envelope of the epidemic. 
The branching envelopes of the simple mean-field epidemics are ordinary Galton- 
Watson processes; the envelopes of the spatial epidemics are branching random 
walks. There is a natural coupling of each epidemic with its branching envelope in 
which the set of infected individuals in the epidemic is at each time (and in the 
spatial models, at each location) dominated by the corresponding set of individuals 
in the branching envelope. 
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Fig 2. 



Following is a detailed description of the natural coupling of the simple SIS 
epidemic with its branching envelope. The branching envelope is a Galton- Watson 
process Z n with offspring distribution Binomial- (N,p), where p is the infection 
parameter of the epidemic, and whose initial generation Zq coincides with the set 
of individuals who are infected at time 0. Particles in the Galton- Watson process 
are marked red or blue: red particles represent infected individuals in the coupled 
epidemic, while blue offspring of red parents represent attempted infections that are 
not allowed because the attempt is made on an individual who is not susceptible, 
or has already been infected by another contagious individual. Colors are assigned 
as follows: (1) Offspring of blue particles are always blue. (2) Each red particle 
reproduces by tossing a p— coin N times, once for each individual i in the population. 
Each Head counts as an offspring, and each represents an attempted infection. If 
several red particles attempt to infect the same individual i, exactly one of these is 
marked as a success (red), and the others are marked as failures (blue). Also, if an 
attempt is made to infect an individual who is not susceptible, the corresponding 
particle is colored blue. Clearly, the collection of all particles (red and blue) evolves 
as a Galton- Watson process, while the collection of red particles evolves as the 
infected set in the SIS epidemic. See the figure below for a typical evolution of the 
coupling in a population of size N — 80, 000 with p = 1/80000 and 200 individuals 
initially infected. 

2. Critical behavior: mean-field case 

When studying the behavior of the simple SIR and SIS epidemics in large pop- 
ulations, it is natural to consider the scaling p — p^ = Xn/N for the infection 
parameter p. In this scaling, A = An is the mean of the offspring distribution in 
the branching envelope. If A < 1 then the epidemic will end quickly, even if a large 
number of individuals are initially infected. On the other hand, if A > 1 then with 
positive probability (approximately one minus the extinction probability for the 
associated branching envelope), even if only one individual is initially infected, the 
epidemic will be large, with a positive limiting fraction of the population eventually 
being infected. The large- N behavior of the size of the SIR epidemic in this case is 
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well understood: see for example [12] and 

2. 1 . Critical scaling: size of the epidemic 

The behavior of both the SIS and SIR epidemics is more interesting in the critical 
case Ajy ~ 1. When the set of individuals initially infected is sufficiently small rel- 
ative to the population size, the epidemic can be expected to evolve in much the 
same manner as a critical Galton- Watson process with Poisson-1 offspring distri- 
bution. However, when the size of the initially infected set passes a certain critical 
threshold, then the epidemic will begin to deviate substantially from the branching 
envelope. For the SIR case, the critical threshold was (implicitly) shown by [U| and 
[H (see also 12 1) to be at N 1 / 3 , and that the critical scaling window is of width 

AT-V3. 

Theorem 1 ([IH, fl). Assume that p N = 1/N + a/TV 4 / 3 + o(n~ 4 / 3 ), and that the 
number J$ of initially infected individuals is such that Jq /iV 1 / 3 — > b as the popu- 
lation size N — > oo. Then as N — ► co, the size Un '■= Ylt ^ t °^ e 2/ s ^ e asymptotic 
law 

(4) U N /N 2 / 3 T h 

where is the first passage time to the level b by Wt + t 2 /2 + at, and Wt is a 
standard Wiener process. 

The distribution of the first passage time Tf, can be given in closed form: See 



111], also |8|, [13|. 

For the critical SIS epidemic, the critical threshold is at N 1 ^ 2 , and the critical 
scaling window is of width N~ 3 / 2 : 

Theorem 2 ([4]). Assume that p N = 1/N + a/N 3 / 2 + o(n~ 3 / 2 ), and that the 
initial number of infected individuals satisfies ~ bN 1 / 2 as N — > co. Then the 
total number of infections Un '■= J3 t Jt during the course of the epidemic obeys 

(5) U N /N r(6 - o; -o) 

where r(x;y) is the time of first passage to y by a standard Ornstein-Uhlenbeck 
process started at x. 



2.2. Critical scaling: time evolution of the epidemic 

For both the SIR and SIS epidemics, if the number of individuals initially infected 
is much below the critical threshold then the evolution of the epidemic will not differ 
noticeably from that of its branching envelope. It was observed by [3] (and proved 
by @) that a (near-) critical Galton- Watson process initiated by a large number 
M of individuals behaves, after appropriate rescaling, approximately as a Feller 
diffusion: In particular, if Z^f is the size of the nth generation of a Galton- Watson 
with Zq ~ bM with offspring distribution Poisson(l + a/M)then asM-^ co, 

(6) Z™ t] /M^Y t 
where Y t satisfies the stochastic differential equation 



(7) 



dY t = aY t dt+ y/YtdWt, 
Yq = b. 
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What happens at the critical threshold, in both the SIR and SIS epidemics, is that 
the deviation from the branching envelope exhibits itself as a size-dependent drift 
in the limiting diffusion: 

Theorem 3 (i). Let J N ( n) = Jr^i be the number infected in the nth generation 
of a simple SIS epidemic in a population of size N . Then under the hypotheses of 
Theorem [H 



(8) J N (VNt)/VN -Z+Y t 
where Yq — b and Y t obeys the stochastic differential equation 

(9) dY t = (aY t - Y t 2 ) dt + y/YtdWt 

Note that the diffusion ([9]) has an entrance boundary at oo, so that it is possible 
to define a version Yt of the process with initial condition Yq — 0. When the SIS 
epidemic is begun with ^> v^V initially infected, the number infected will 
rapidly drop (over the first e-\/~N generations) until reaching a level of order y/~N, and 
then evolve as predicted by {5| . The following figure depicts a typical evolution in a 
population of size iV = 80, 000, with infection parameter p = 1/N and Iq = 10, 000 
initially infected. 

Theorem 4 ([4]). Let J N (n) = JP, and R N (n) = be the numbers of infected 
and recovered individuals in the nth generation of a simple SIR epidemic in a 
population of size N. Then under the hypotheses of Theorem^ 

(N-V*J N (NVH)\ [J(t)\ 

[W) {n- 2 / 3 r n (n 1 /H)J \R(ty 

where J — b, R a — 0, and 

(11) dJ(t) = (aJ(t) - J{t)R(t)) dt + y/J(t)dWt, 

dR(t) = J(t) dt. 

Theorems HH2] can be deduced from Theorems |3H1] by simple time-change argu- 
ments (see |4j). 



Fig 3. 
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The critical thresholds for the SIS—d and SIR—d epidemics can be guessed by sim- 
ple comparison arguments using the standard couplings of the epidemics with their 
branching envelopes. Consider first the critical SIS epidemic in a population of size 
N. Recall (Section II .4[) that the branching envelope is a critical Galton- Watson 
process whose offspring distribution is Binomial- (N, 1/N). The particles of this 
Galton- Watson process are marked red or blue, in such a way that in each genera- 
tion the number of red particles coincides with the number of infected individuals 
in the SIS epidemic. Offspring of blue particles are always blue, but offspring of 
red particles may be either red or blue; the blue offspring of red parents in each 
generation represent attempted infections that are suppressed. 

Assume that initially there are N a infected individuals and thus also N a indi- 
viduals in the zeroth generation of the branching envelope. By Feller's theorem, we 
may expect that the extinction time of the branching envelope will be on the order 
N a , and that in each generation up to (shortly before) extinction the branching 
process will have order N a individuals. If a is small enough that the SIS epidemic 
obeys the same rough asymptotics (that is, stays alive for 0(N a ) generations and 
has 0(N a ) infected individuals in each generation), then the number of blue off- 
spring of red parents in each generation will be on the order N x (N 2a /N 2 ) (because 
for each of the N individuals of the population, the chance of a double infection 
is about N 2a /N 2 ). Since the duration of the epidemic will be of the same rough 
order of magnitude as the size of the infected set in each generation, there should 
be at most 0(1) blue offspring of red parents in any generation (if there were more, 
the red population would die out long before the blue). Thus, the critical threshold 
must be at A 1 / 2 . 

A similar argument applies for the SIR epidemic. The branching envelope of the 
critical SIR is once again a critical Galton- Watson process with offspring distribu- 
tion Binomial- (N, 1/N), with constituent particles again labeled red or blue, red 
particles representing infected individuals in the epidemic. The rule by which red 
particles reproduce is as follows: Each red particle tosses a p— coin N times once 
for each individual i in the population. Each Head counts as an offspring, and rep- 
resents an attempted infection. However, if a Head occurs on a toss at individual 
i where i was infected in an earlier generation, then the Head results in a blue 
offspring. Similarly, if more than one red particle tosses a Head at an individual i 
which has not been infected earlier, then one of these is labeled red and the excess 
are all labeled blue. 

Assume that initially there are A a infected individuals. As before, we may expect 
that the extinction time of the branching envelope will be on the order N a , and 
that in each generation up to extinction the branching process will have order 
N a individuals. If a is small enough, the extinction time and the size of the red 
population will also be 0(N a ). Consequently, the size of the recovered population 
will be (for all but the first few generations) on order N 2a . Thus, in each generation, 
the number of blue offspring of red parents will be on order (N 2a /N) x N a (because 
the chance that a recovered individual is chosen for attempted infection by an 
infected individual is 0(N a /N)). Therefore, by similar reasoning as in the SIS 
case, the critical threshold is at A 1 / 3 , as this is where the the number of blue 
offspring of red parents in each generation is O(l). 
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3. Critical behavior: SIS-1 and SIR-1 Spatial epidemics 

Consider now the spatial SIS-d and SIR-d epidemic models on the d-dimensional 
integer lattice Z d . Assume that the village size N x = N is the same for all sites x £ 
Z d , and that the infection probabilities a(x,y) are nearest neighbor homogeneous, 
and uniform, that is, 



(12) a(x,y) 




if \x-y\ <1; 
otherwise. 



3.1. Scaling limits of branching random walks 

The branching envelope of a spatial SIS—d or SIR—d epidemic is a nearest neighbor 
branching random walk on the integer lattice Z d . This evolves as follows: Any parti- 
cle located at site x at time t lives for one unit of time and then reproduces, placing 
random numbers £ y of offspring at the sites y such that \y — x\ < 1. The random 
variables £ y are mutually independent, with Binomial- (N,pn) distributions. 

The analogue for branching random walks of Feller's theorem for Galton- Watson 
processes is Watanabe's theorem. This asserts that, after suitable rescaling, as the 
particle density increases, critical branching random walks converge to a limit, 
the Dawson- Watanabe process, also known as super Brownian motion. A precise 
statement follows: Consider a sequence of branching random walks, indexed by 
M = 1, 2, . . ., with offspring distribution Binomial- (N,pm) as above, and 

la 

(13) Pm=Pn,m 



(2d+l)N NM 

(Note: N may depend on M.) The rescaled measure-valued process X^ 1 associated 
with the Mth branching random walk puts mass 1/M at location xj\[M and time t 
for each particle of the branching random walk that is located at site x at time [Mi] . 
(Note: The branching random walk is a discrete-time process, but the associated 
measure- valued process runs in continuous time.) 

Watanabe's theorem ((HI). Assume that the initial values Xq 1 converge weakly 
( as finite Borel measures on W 1 ) to a limit measure Xq . Then under the hypothesis 
(|13|) the measure-valued processes X t M converge in law as M — > oo to a limit X t : 

(14) X t M =>X t . 

The limit process is the Dawson- Watanabe process with killing rate a and ini- 
tial value Xq. (The term killing rate is used because the process can be obtained 
from the "standard" Dawson- Watanabe process (a = 0) by elimination of mass 
at constant rate a.) The Dawson- Watanabe process X t with killing rate a can be 
characterized by a martingale property: For each test function <fi G C^(R d ), 

(15) {X t , <p) - (A , 0) - - / (X„A4>) ds + a (X s , if) ds 

* Jo Jo 

is a martingale. Here a 2 = 2d/(2d + 1) is the variance of the random walk ker- 
nel naturally associated with the branching random walks. It is known fl~oj ] that 
in dimension d = 1 the random measure X t is for each t absolutely continuous 
relative to Lebesgue measure, and the Radon-Nikodym derivative X(t, x) is jointly 
continuous in t, x (for t > 0). In dimensions d>2 the measure X t is almost surely 
singular, and is supported by a Borel set of Hausdorff dimension 2 
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3.2. Spatial SIS-1 and SIR-1 epidemics: critical scaling 

As in the mean- field case, there are critical thresholds for the SIS-1 and SIR-1 
epidemics at which they begin to deviate noticeably from their branching envelopes. 
These are at N 2 / 3 and iV 2//5 , respectively: 

Theorem 5 (|5[). Fix a > 0, and let X^ be the discrete-time measure-valued 
process obtained from an SIS— 1 or an SIR— 1 epidemic on a one- dimensional grid 
of size- N villages by attaching mass 1/N a to the point (t,x/N a / 2 ) for each infected 
individual at site x at time [tN a ]. Assume that Xff converges weakly to a limit 
measure Xq as the village size N — > oo. Then as N — > oo, 

(16) X { N Nat] ^X t , 

where X t is a measure- valued process with initial value Xq whose law depends on 
the value of a and the type of epidemic ( SIS or SIR ) as follows: 

(a) SIS: Ifa<l then X t is a Dawson- Watanabe process with variance a . 

(b) SIS: If a = | then X t is a Dawson- Watanabe process with variance a 2 and 
killing rate 

(17) 9(x,t) = X(x,t)/2. 

(c) SIR: // a < | then X t is a Dawson- Watanabe process with variance a 2 . 

(d) SIR: // a = | then X t is a Dawson- Watanabe process with variance a 2 and 
killing rate 



(18) 



6{x, t) = X(x, t) [ X(x,s)ds. 
Jo 



The Dawson- Watanabe process with variance a 2 and (continuous, adapted) 
killing rate 0(t,x,uj) is characterized [2] by a martingale problem similar to (| 1 5|) 
above: for each test function tfi £ C 2 (M.), 

(19) {X t ,4>) - {X ,4>) - - / (X Sl A<f>)ds+ (X s ,9<p)ds 

z Jo Jo 

is a martingale. The law of this process is mutually absolutely continuous relative 
to that of the Dawson- Watanabe process with no killing, and there is an explicit 
formula for the Radon-Nikodym derivative - see 



3.3. Critical scaling for spatial epidemics: heuristics 

Arguments similar to those given above for the mean-field SIS and SIR epidemics 
can be used to guess the critical thresholds for the spatial SIS-1 and SIR-1 epi- 
demics. For the spatial epidemics, the associated branching envelopes are branching 
random walks. In the standard couplings, particles of the branching envelope are 
labeled either red or blue, with the red particles representing infected individuals 
in the epidemics. As in the mean-field cases, offspring of blue particles are always 
blue, but offspring of red particles may be either red or blue; blue offspring of red 
parents represent attempted infections that are suppressed. These may be viewed 
as an attrition of the red population (since blue particles created by red parents 
are potential red offspring that are not realized!). 
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Consider first the SIS-1 epidemic. Assume that initially there are N a parti- 
cles, distributed (say) uniformly among the TV"/ 2 sites nearest the origin. Then by 
Feller's limit theorem (recall that the total population size in a branching random 
walk is a Galton- Watson process), the branching envelope can be expected to sur- 
vive for Op(N a ) generations, and at any time prior to extinction the population 
will have Op(N a ) members. These will be distributed among the sites at distance 
Op{N a ' 2 ) from the origin, and therefore in dimension d = 1 there should be about 
Op(N a / 2 ) particles per site. Consequently, for the SIS—1 epidemic, the rate of 
attrition per site per generation should be Op(N a ~ 1 ), and so the total attrition 
rate per generation should be Op(N 3a ^ 2 ^ 1 ). If a — 2/3, then the total attrition 
rate per generation will be Op(l), just enough so that the total attrition through 
the duration of the branching random walk envelope will be on the same order of 
magnitude as the population size N a . 

For the SIR—1 epidemic there is a similar heuristic calculation. As for the SIS—1 
epidemic, the branching envelope will survive for Op(N a ) generations, and up 
to the time of extinction the population should have Op(N a ) individuals, about 
Op(N a / 2 ) per site. Therefore, through N a generations, about N a x N a / 2 numbers j 
will be retired, and so the attrition rate per site per generation should be Op(N a / 2 x 
N 3a / 2 ), making the total attrition rate per generation Op(N 5a / 2 ). Hence, if a — 2/5 
then the total attrition per generation should be Op(l), just enough so that the 
total attrition through the duration of the branching random walk envelope will be 
on the same order of magnitude as the population size. 

3-4- Critical scaling in dimensions d > 2 

In higher dimensions, the critical behavior of the SIS-d and SIR-d epidemics appears 
to be considerably different. We expect that there will be no analogous threshold 
effect, in particular, we expect that the epidemic will behave in the same manner as 
the branching envelope up to the point where the infected set is a positive fraction of 
the total population. This is because in dimensions d>2, the particles of a critical 
branching random walk quickly spread out, so that (after a short initial period) 
there are only Op(l) (in dimensions d > 3) or Op (log iV) (in dimension d = 2) 
particles per site. (With N particles initially, a critical branching random walk 
typically lives 0(N) generations, and particles are distributed among the sites at 
distance 0{s/N) from the origin; in dimensions d > 2, there are 0(N d / 2 ) such sites, 
enough to accomodate the O(N) particles of the branching random walk without 
crowding.) Consequently, the rate at which "multiple" infections are attempted 
(that is, attempts by more than one contagious individual to simultaneously infect 
the same susceptible) is only of order Op(l/N) (or, in dimension d = 2, order 
OpQogN/N)). 

The interesting questions regarding the evolution of critical epidemics in dimen- 
sions d > 2 center on the initial stages, in the relatively small amount of time (order 
o(N) generations) in which the particles spread out from their initial sites. These 
will be discussed in the forthcoming University of Chicago Ph. D. dissertation of 
Xinghua Zheng. 

3.5. Weak convergence of densities 

There is an obvious gap in the heuristic argument of Section [3.31 above: Even if the 
total number of infected individuals is, as expected, on the order N a , and even if 
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these are concentrated in the sites at distance on the order N a / 2 from the origin, it 
is by no means obvious that these will distribute themselves uniformly (or at least 
locally uniformly) among these sites. The key step in filling this gap in the argument 
is to show that the particles of the branching envelope distribute themselves more 
or less uniformly on scales smaller than N a / 2 . 

Consider, as in Section 13.11 a sequence of branching random walks, indexed by 
M = 1, 2, ... , with offspring distribution Binomial- (N,pm) as above, and pm given 
by (fl3|) . Let Y t M (x) be the number of particles at site x at time [t], and let X M (t, x) 
be the continuous function of t > and x € R obtained by linear interpolation from 
the values 



(20) x M (t, x) = Y Mt{y/Mx) for m g and ^- Ix g ^ 



Theorem 6 (@j). A ssume that d = 1. Assume also that the initial particle config- 
uration is such that all particles are located in an interval [—KyM, ka/M] and such 
that the initial particle density satisfies 

(21) X M {0,x) =*> X{0,x) 

as M — > oo for some continuous function X(0,x) with support [—k,k]. Then as 
M -> oo, 

(22) X M (t,x) =► X{t,x), 

where X(t,x) is the density function of the Dawson- Watanabe process with killing 
rate a and initial value X(0,x). The convergence is relative to the topology of uni- 
form convergence on compacts in the space C(R+ x R) of continuous functions. 

Since the measure-valued processes associated with the densities X M (t,x) are 
known to converge to the Dawson- Watanabe process, by Watanabe's theorem, to 
prove Theorem [6] it suffices to establish tightness. This is done by a somewhat 
technical application of the Kolmogorov-Chentsov tightness criterion, based on a 
careful estimation of moments. See [5j for details. 

It is also possible to show that convergence of the particle density processes holds 
in Theorem [5] 
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